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ION  TEMPERATURE  SENSITIVE  EFFECT  IN 
TRANSIENT  LANGMUIR  PROBE  RESPONSE 
by 

Juan  R.  Sanmartln 

Department  of  Mechanical  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 

ABSTRACT 

A  theory  is  presented  for  a  method,  recently  proposed  by  Hester 
and  Sonin,  of  determining  the  ion  temperature  in  a  plasma  by  measur¬ 
ing  the  transient  current  response  of  a  highly  negative,  cylindrical 
Langmuir  probe  under  conditions  where  collection  is  collision-free 
and  the  ratio  of  probe  radius  to  Debye  length  is  small.  The  ion 
component  of  the  current  does  not  approach  its  final,  steady-state 
value  monotonlcally,  but  exhibits  a  strong,  ion  temperature  dependent 
overshoot  in  the  first  few  ion-plasma  periods  following  the  biasing 
of  the  probe.  Analytical  formulae  are  derived  for  a  Maxwellian 
plasma,  and  convenient  graphical  results  are  presented.  The  possible 
masking  of  the  overshoot  by  a  transient  displacement  current  is 
discussed  and  found  to  be  often  negligible.  For  other  conditions  an 
alternative  measuring  method  is  suggested  in  which  the  contribution 
of  the  displacement  current  cancels  out. 


I  -  Introduction 


Hester  and  Sonin*  have  recently  proposed  a  method  of  determining 
the  ion  temperature  in  a  plasma  by  measuring  the  transient  current 
drawn  by  a  long,  cylindrical  Langmuir  probe  to  which  a  large,  negative 
potential  has  been  applied  step-wise.  Thev  performed  some  computa¬ 
tions  showing  that  the  ion  current  exhibits  an  initial  overshoot  which 
is  very  sensitive  to  the  ion  temperature  when  the  probe  radius  is 
small  compared  with  the  Debye  length.  Such  an  effect  would  be  of 
considerable  Interest  for  diagnostic  purposes  because  steady-state 

probe  response  is,  under  all  conditions,  very  weakly  dependent  on  that 
2 

plasma  property  . 

In  this  paper,  we  present  a  theoretical  formulation  of  that 
problem  which  yields  both  a  simple  analytical  formula  for  the  tran¬ 
sient  current  and  a  clear  understanding  of  the  overshoot.  The  treat¬ 
ment  is  limited  to  the  first  ion-plasma  period  following  the  biasing 
of  the  probe,  which  is  the  region  where  the  response  is  most  sensitive 
to  the  ion  temperature.  The  basic  approach  is  discussed  in  the  next 
section.  In  Sec.  Ill  we  derive  some  needed  results  about  the  potential 
field  set  up  by  the  probe.  The  phenomenon  of  an  ion  current  overshoot 
is  discussed  in  Sec.  IV  and  in  Sec.  V  a  formula  is  derived  for  the 
current  and  convenient  graphical  results  are  presented.  Finally,  a 
discussion  of  experimental  aspects  of  the  method,  and  its  usefulness, 
is  given  in  Sec.  VI.  In  particular,  the  possible  masking  of  the  over¬ 
shoot  by  the  displacement  current,  which  is  registered  in  the  external 
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measuring  circuitry  together  with  the  ion  current,  is  discussed.  It 
is  shown  that  such  masking  is  often  negligible.  For  other  conditions 
where  it  is  not  negligible,  an  alternative, pulsed  probe  method  is 
suggested  in  which  displacement  effects  cancel  out. 

II  -  Basic  formulation 

Consider  the  problem  of  a  quiescent  plasma  and  an  infinitely 
long,  cylindrical  Langmuir  probe  of  radius  r^  to  which  a  negative 
potential  (relative  to  the  plasma)  is  applied  at  t  ■  0.  It  is 
desired  to  obtain  the  ion  current  to  the  probe  for  positive  times. 

The  potential  field,  V,  is  governed  by  Poisson's  equation,  which  in 
nondimens ional  form  reads 


13  3  . 

p  3p  P  3p  ^ 


2 

c 


with  boundary  conditions 


(1) 


<Ji  (p  *  1)  *  <{»  =  -  eV  /kT 

P  P  e 

“  0  , 

where 


0  *  -  eV/<T  ,  o  ■  r/r 
e  p 

e  -  rp/*D  »  n*  *  ZiNi/No 


n 

e 


Vo 
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and  are  the  electron  and  Ion  densities,  Tg  and  Nq  the  unperturbed 

electron  temperacure  and  density,  r  is  the  radial  distance  to  the  probe 

axis,  the  ion  charge  number  and  X^  the  electron  Debye  length, 

2  1/2 

(>cT  /4ire  N^)  .  It  will  be  here  assumed  that  all  mean  free  paths  are 

large  compared  with  X^,  so  that  if  c  $  0(1),  current  collection  will  be 
collision-free  and  ng  and  n^  will  be  given  by  the  corresponding  time- 
dependent  Vlasov  equation.  Simple  dimensional  considerations  then 
show  that 


-1(T  >  0) 

J. 


e,  *p. 


S, 


(2) 


where 


W. 


T  ■  U  t 

P 


j  is  the  ion  current  density,  j  the  steady-state  current  density 

corresponding  to  the  probe  potential  V  ,  u>p  the  ion-plasma  frequency 
2  1/2 

(4xe  Z.N_/m.)  ,  m  the  electron  mass,  and  m,  and  T.  the  ion  mass 

i  0  i  e  11 

and  temperature  respectively;  Eq.  (2)  will  also  depend  on  the  probe 
potential  for  negative  times,  if  that  is  not  zero.  One  expects 
j/jm  +  1  as  i  +  •. 

The  general  solution  of  the  problem  defined  above  is  quite 
difficult.  We  are  interested,  however,  in  specific  conditions  under 
which  j  /.I  does  not  increase  monotonically  from  j(x  ■  0)/.)^  to  1,  but 
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exhibits  a  strong  overshoot.  It  Is  quite  fortunate  that  such  conditions 
result  also  in  a  considerable  simplification  of  the  problem.  We  assume 
the  following: 

(a)  The  probe  is  at  plasma  potential  at  t  <  0.  The  selec¬ 
tion  of  any  other  initial  probe  potential  would  somewhat  mask  the  in¬ 
fluence  of  on  the  overshoot,  since  this  influence  depends  critically 
on  the  initial  distribution  function  in  the  neighborhood  of  the  probe. 

However,  the  above  condition  makes  the  analysis  easier  because  of  the 

2 

simplicity  of  the  initial  plasma  state;  we  have  in  particular 

<Kp,  t  <  0)  -  0  , 

n  (p,  t  $  0)  ■  n.  (p,  t  $  0)  ■  1  -  it  1  sin  1  p  1  ,  (3) 

e  l 

](i  f  0)  ■  e  N0(2tt)"1/2  (KTi/mj)in  .  (4) 

(b)  ijip  >>  1.  This  assumption  allows  us  to  neglect  the  elec¬ 
tron  current.  Moreover,  it  leads  to  a  simple  expression  for  ng, 

ng  -  exp  (-  il<)  (5) 

valid  everywhere  except  in  the  neighborhood  of  the  probe  where  ng  will 
be  exponentially  small  anyway  [see  also  condition  (d)  below].  The 


equilibrium  with  the  field  will  be  reached  in  a  few  electron-plasma 
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periods,  so  chat  the  dependence  of  (2)  on  m  / m.  will  disappear  for,  say, 

e  l 

t  >  0.1. 

(c)  ip  >>  8,  and  (d)  e  <<  1.  These  two  conditions  are 
P 

found  to  be  essential  for  the  appearance  of  a  large  overshoot.  Con¬ 
dition  (d)  also  implies  that  steady-state  current  collection  will  be 

2 

orbital  motion  limited,  so  that 

-1»  “  e  N0  2l/2'~1  ZieVp/mi)1/2  ’  (6) 

2 

and,  assuming  the  unperturbed  plasma  to  he  Maxwellian, 

1  f  .1  l  'lln  ^  P®  l/o 

n^i  ■*  -)  -  1  -  ±  j  ds  exp(-  s)  sin  1  [±  (-*  y  p  '  1  • (7) 

0 

Furthermore,  (d)  greatly  simplifies  the  determination  of  ♦  ,  and  con¬ 
tributes  also  to  the  validity  of  Eq.  (5),  since  a  region  of  order  of 

rp,  where  Eq.  (5)  falls,  then  contains  a  negligible  amount  of  negative 
charge  even  if  ng  -  0(1). 

All  these  conditions  greatly  simplify  the  analysis.  First,  it 

is  possible  to  derive  an  accurate  analytical  expression  for  3iJ//3p 

without  solving  simultaneously  the  ion  Vlasov  equation.  That 

expression  is  only  valid  in  a  certain  (large)  neighborhood  of  the  probe. 

This  restricts  the  theory  to  times  bound  by  an  upper  value  t  ,  since 

m 

for  sufficiently  long  times  the  probe  will  be  collecting  ions  that 
were  initially  outside  that  neighborhood.  Fortunately,  again,  the  over¬ 
shoot  occurs  at  shorter  times,  and  is,  in  fact,  most  sensitive  to  T^  in 
the  T-range  where  the  theory  is  valid. 
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Second,  once  the  field  Is  known  ion  trajectories  may  be  determined 
explicitly.  Moreover  the  current  to  the  probe  is  now  linear  in  the 
initial  ion  distribution  function  and  this  allows  a  further  simplifica¬ 
tion.  Let  f^  be  the  ion  distribution  function  after  integration  over 
the  velocities  parallel  to  the  probe.  We  shall  assume  that  far  away 
from  the  probe  f^  is  isotropic 

^(p  -*■“)■  Z"1  Nq  g(vA) 

where  vA  Is  the  ion  speed  in  the  p-$  plane  ($  being  a  polar  angle) . 
f^  is  normalized  such  that 

90 

|  2nvA  g(vA)  dvA  -  1 
0 

for  a  Maxwellian  plasma 

2 

m^  -  nijVA 

8  "  2wieT^  exp  (  2kT1 
Then  we  clearly  have 

f^x  $  0)  -  Z~\  g(vA) 

-  0  , 


(8) 

,  |y|  >  sin  1  p  1  (9a) 

M  <  sin"1  p”1  (9b) 


| 
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where  y  is  a  velocity-vector  angle  (-  n  s  y  $  it)',  tills  is  illustrated 
In  Fig.  1.  For  any  g,  Eqs.  (9)  result  in  the  initial  density  given 
by  Eq.  (3);  the  ion  depletion  is  due  to  the  particle  loos  caused  by 
thermal  motion.  Now  let  j*  be  the  current  density  obtained  by  assum¬ 
ing  an  initial  ion  distribution  function  with  all  ions  having 
velocities  with  the  same  direction  and  magnitude,  v*;  see  wig.  2, 
where  the  polar  axis  has  been  chosen  parallel  to  the  velocity  for 
simplicity.  (We  note  that  in  the  shaded  area  there  are  no  ions  while 
outside  it  the  density  is  uniform,  N^/Z^.)  Then  one  easily  concludes 
from  the  symmetry  of  the  problem  that  j  may  be  written  as 

00 

j(i  *  0)  »  J  2x  v*g(vA*)  dv*  j*(v*,  T  i  0)  .  (10) 

0 

Finally  the  determination  of  j*  itself  may  be  simplified  by 
using  the  following  approach.  The  non-shaded  region  of  Fig.  2  may 
be  divided  in  two  mutually  exclusive  parts.  A*  and  B*.  A  point 
o— «t>  belongs  to  A*  if  an  ion  having  such  initial  coordinates  will  hit 
the  probe  at  some  time  t  >  0.  Computing  the  time  of  flight  of  an 
A*-ion  to  the  probe  we  may  then  define  a  function  a*(r)  that  repre¬ 
sents  the  area  of  region  A*  that  has  been  "collected"  by  the  time  t. 
From  the  uniformity  of  the  density  it  follows  that 

2irr  j*  •  eN-  da*/dr 
P  0 


or 
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J*  .  e  dh* 

(8i|/  )1/2  dT 

P 

.  »  2 
where  a*  is  a  nondimensional  area,  a*  -  a*/r  .  For  the  time-averaged 

P 

current  we  have 


i *  =  1 


0 


dr' 


ta* 


«v 


1/2 


(ID 


In  the  next  section  we  will  find  the  aforementioned  expression 
for  the  electric  field,  3i|>/3p.  In  See.  IV  we  will  determine  the 
boundaries  of  the  region  A*  and  discuss  the  overshoot  phenomenon. 

In  Sec.  V  we  will  compute  the  area  a*  and  from  there  the  current  to 
the  probe. 

Ill  -  The  potential  field 

In  the  next  two  sections  we  shall  have  to  integrate  the  electric 
field  3iJj/3p  along  the  trajectory  of  a  particle  to  obtain  its  velocity 
and  its  time  of  flight  to  the  probe.  A  formal  expression  for  3<|>/3p 
can  be  obtained  from  Eq.  (1), 


3jj» 

3p 


-  ♦_  S(t)  .2  ? 

- E_ - +  1  p'  dp'  v (p  ' ,  t) 

P  P  J 

1 


(12) 


where 


6 


3ip/ 3p  |  , 

|  P-1 


* 


Eq.  (12)  may  be  rewritten  as 
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3jJ/ 
3  P 


*n  <5(0  2  2, 

- E-- - +  —  - <  v  (p  ,  T  )  > 

p  p  l 


(13) 


<v>  being  given  by 

(o2  -  l)<v> 


v  d(p'2) 


(14) 


We  note  Chat  6  is  proportional  to  the  electric  field  at  the  probe. 
Its  relation  to  the  overall  charge  density  follows  easily  from  (12): 


.-1  2 


w  w 

j  p  dp  v  ■  (j 


P  dp  v)/( 


p  dp  v  In  p) 


(15) 


Thus  6  1  represents  the  logarithm  of  a  characteristic  shielding 
distance. 

To  carry  out  the  aforementioned  Integrations  we  shall  substitute 
constants  6  and  v  for  6(t)  and  <v(p,  t)>  in  Eq.  (13).  Our  purpose 
now  is  to  show  that  <5  and  (with  some  restrictions  on  p)  <v>  vary 
little,  and  to  find  close  bounds  for  them;  this  will  allow  us  to  choose 
appropriate  values  for  6  and  v,  and  to  determine  the  error  in  the 
approximations  used.  The  basic  approach  will  be  to  study  the  limits 
t  ■*  0  and  t  ■*  °»,  for  which  v  is  known  explicitly  (see  Eqs.  (3),  (5) 
and  (7)J. 

3 

(A)  Bounds  for  6  at  t  Jfc  0  and  t  -*>  “. 

We  must  first  establish  the  following  theorem:  Let  tl^vo) 


be  defined  by 
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If  0  “  Fk^k(p)* 01  031  *  (16) 

<»k(l)  ■  ^  .  Ij»k(")  *  0  (17) 

where  k  »  1,  2  and  is  bounded,  compatible  with  Eqs.  (17),  and 
such  that  ^k  is  continuous  and  monotonlc.  Let 

F2(f,  ol  5  F1[f,  p]  Oil  ,  ,  (18) 

and  F^f,  p]  be  nondecreasing  in  f.  Then,  if  6^  =  -  <J>  1  3iJ>k/3p  ■  , 

P  ( o—l 

62  31  61 

To  prove  this  let  us  assume  that  62  <  6^  Integrating  (16)  twice 
we  get 

o 

tf'k  •  •I'pll  -  6k  In  p]  +  |  o’  dp'  In  F[i|<k  (o'),  p’J 

1 

Thus  in  a  certain  neighborhood  of  1,  we  have  i^2  >  Let  s  be  the 

minimum  value  of  p  such  that  i|»2(s)  ■  'l'^(s);  such  value  exists,  whether 
finite  or  infinite,  because  ti^O*)  •  <J^(“).  It  follows  that 
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8 

0  •  i}>^(s)  -  ij>2(s)  "  *n  8  +  j  P  dp  In  -* 

1 

*  ^[^(p),  pj  -  F2U2(p),  p )  >  $  ^(6^  -  62)  In  8 
8 

+  |  P  dp  In  5-  (F1[*1(p),  pi  -  F1[<,2(p),  p]>  <  0  , 

1 

proving  Chat  62  cannot  be  less  than  6^. 

Upper  and  lower  bounds  for  6(t  %  0)  can  now  be  found.  The  right- 
hand  side  of  Eq.  (1)  is  then,  according  to  Eqs.  (3)  and  (5), 


F  “  e2[l  -  exp  (-  \p)  3  ~  ~  sin  *  — 

*  p 


Introducing 


F^ul^p),  p]  -  e2[l  -  exp(-  <|>  )J  -  ~  sin”1  ^  j  1  -  exp(-  *  )  , 

L  L  p  X  P  2  P 


2,(1)  c  ,  -1  1 
c  *2  -  ~  ®in  - 


Si-  exp(-  i/i  )  , 

t  P 


our  theorem  can  be  seen  to  apply  if  F  -*•  F^,  if>  -*•  and  F^1^  F2 
■*  i >2«  Therefore 


6<t  %  0)  $  6<1} 


(19) 


Since  F 


(1) 


is  linear,  5^  is  easily  found;  neglecting  exp(-  ^p), 
and  smaller  terms  in  favor  of  unity,  we  get 


I 


-  [1  -  ((T1  (1  -  AqjM)]  In  o (1)  . 


(20a) 


is  such  that  ^  '[oqj]  ■  1  -  exp(-  and  is  given  by 


*  . ,  e2pq)  +. f-e2oli)  +  CQ  ki(cp(d) 
*P  1  A  +  1  2  £0(1)  K0(cp(1)> 


n  {P(1)  + 


I  ^ 

Ko 


K.(cp)  dp 


(ep(1)) }1  ln  p(l) 


(20h) 


For  a  lower  bound,  let  us  introduce 


lfV2)(p>,  p]  ■  c 2 [  1  -  exp(-  p)]  -  sin 


e  ,-11  (2) 

-  — -  sin  *.  i  P 

IT  pi 


:2  AO  .  c_  8ln-l  I 

P  1  IT  P 


*1  S  P  , 


where  we  wrote  p  for  1  -  exp(-  to  simplify  the  equations.  Our 

(2)  (2) 

theorem  applies  again  with  F  ■*  F^t  &  -*■  i|*2  and  -►  F^,  •+ 


6(x  %  0)  $  6. 


(2) 

5^  is  given  by 


l/6^2)  -  [1  -  Vp1  (1  -  e2o22)/4)]  [ln  p(2)  +  0.23]  ,  (22a) 
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,  e2p(2)  fe2p(22  K1(CP(2)) 

VP  1  ‘  4  1  2  +  Cp (2)  K0(cp(2)T 


“  r  1‘26  <p(2)  + 


r  K0(ep)  dp 

J  K^rmin  °(2) +  °-23'  • 


(2) 


For  6(t  ■  •),  we  have  [according  to  (7)] 


7  f  .  *  +  8s 

F  -  e  [1  -  exp(-  i|»)]  -  —  J  ds  exp(-  s)  sin  l“(^jTg“*)  ' 


as  the  right-hand  side  of  Eq.  (1).  An  upper  hound  follows  imnediately 


6(t  -  •»)  S  6(t  %  0)  $  6^ 


because  3F/3i l>  >  0  in  (23)  and 


1  (  ,  ,  .  .  -1,1  ,*p  +  Ss,l/2,  .1  ,-11 

V  1  ds  exp(-  5)  ,ln  I'o'  <* +-BV')  1  3  T  sin  0 


for  0  $  \h  i  •b  .  For  a  lower  bound,  we  note  that  for  6  >  8 
P 


6(t  -  «,  8)  >.  6(t  -  »,  B) 


since 


iji  +  8s  +  Bs 

.  P____  *  J L__ 

<i>  +  8s 


(22b) 


]  (23) 


(24) 


4>  +  fls 
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for  s  >  0  and  0  $  ip  S  ip  ;  thus,  it  suffices  to  find  a  lower  bound  for 

P 

6 ( t  ■  »,  8*0).  Let  us  Introduce 


-  e2 [  1  -  exp(-  p)]  -  | 


ds  exp(-  8)  [tt,  sin”1  {- 
*  p 


rV:  Savl/2 


P  +  Ss 


> "■}) 


<,<3>  >  P 


(25a) 


-  e2  -1  -  p) 

P 


i*-£l 


ds  exp(-  s)  [x,  sin  1  — 
t  P 


ip  +  8s 
(— * 


(pPp^/p2)  +  8s 


)1/2» 


43)  s  P  , 


(25b) 


where  again  p  •  1  -  exp(-  ijip),  and  in  the  brackets  inside  the  integrals 

above,  the  second  terms  are  to  be  used  whenever  the  arguments  of  sin  1 

are  less  than  1;  otherwise  n/2  should  be  used.  Notice  that  we  have 

(3) 

retained  a  finite  8:  if  we  take  8  0  in  (25),  cannot  be  made  to 

vanish  at  infinity.  However  if  we  first  retain  8  j*  0,  solve  for 

i|»^3^  and  6^3\  and  then  let  8  0,  we  find  a  finite  6^;  the  singularity 

of  this  limit  is  related  to  well  known  differences  between  the  steady- 

state  currents  to  a  negative  cylindrical  probe  under  condition  ■  0 

3  (3) 

and  condition  0.  We  also  note  that  is  not  a  point  function 

but  a  functional  of  because  of  the  appearance  in  Eq.  (25)  of  p 

(3)  (3) 

the  point  where  ■  p;  making  F  ■+•  ip  ^  and  F^  F^, 

(3)  4 

-*■  ip^,  one  can  see  that  (18)  is  nevertheless  satisfied  and  that 
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(3) 

Fj ,  i.e.,  ,  satisfies  a  certain  nondecreasing  condition:  if 

^3^(p')  $  i l>(o')  in  the  entire  range  1  $  p'  $  p,  p]  $  f|3^  [iJ»,o]. 

Accounting  for  this,  our  theorem  may  be  seen  to  apply  again,  giving 


6(t  ■  •,  B)  i  6(t 


6  *  0)  i  { 


(3) 

1 


(26) 


1  / 6 5 3>  »  (1  -  1  -  e2R2/4  +  £2Ri|)1^2  1.26/h}) 

ip  p 


(In  R  +  0.231 


(27a) 


,  e2R2  ,  ,eV  ,  „  Kl(cR>  2  1/2  1.26 

"p  ■  1  -  —  +  >—  +  *  Kjt «  -  E  *p  — 


K.  ( cR)  9  1/?  l  ?<> 

•  {R  +  7k  (dT)  ^  (ln  R  +  °*23)  +  c  V  R  *  *  (27b) 


2 

where  R  ■  P^j/1.26.  We  have  neglected  exp(-  and  c  38 

2 

in  Eqs.  (20)  and  (22),  and  also  <1^/30^. 

In  the  parameter  range  later  considered, 

10“2  10"1  ,  10  $  *  $  102  ,  (28) 

P 


the  terms  neglected  in  deriving  Eqs.  (20),  (22)  and  (27)  result  in  an 

error  of  less  than  1Z.  Within  such  an  error  we  also  find 

l/fij2^  -  1/6^  +  0.23,  and,  within  an  error  of  less  than  1%  for  e  •  10 

-1  (3)  (2) 

and  less  than  2.5Z  for  e  ■  10  ,  we  find  5^  ■  5^  . 
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(B)  Bounds  for  <v>  at  T  %  0  and  T  ■*  «. 
For  T  %  0  we  have 


v 


1 


(-  *) 


(29) 


Since  i|)  >>  1,  we  get  v(l)  %  1/2.  As  p  Increases  v  approaches  unity 

P 

and  remains  close  to  it  until  becomes  small  enough  for  the  last 
term  in  Eq.  (29)  to  be  important.  At  p  *  1.05, v  >  0.60.  At  p  ■ 
as  given  by  Eq.  (20b),  i|»  %  1  since  at  such  distances  thus 

v[p(i)]  %  0.63.  Therefore  v  >  0.60  between  p  ■  1.05  and  p  slightly 
larger  than  p^^;  over  most  of  this  range  v  will  be  close  to  1.  Now, 
substitute  a  constant  6  for  5(t)  in  Eq.  (14)  [see  Eq.  (35)  below], 
and  define 


p 


m 


(2*p  «) 


1/2 


-1 

e 


(30) 


We  have 


3j|> 

3p 


-  \l>  6 

_ JL— 

p 


[1 


2 

For  t|ip  >>  1,  e  <<  1,  Pm  is  quite  large  so  that  the  last  term  in  the 
bracket  can  be  neglected.  The  obvious  requirement  3iJ>/3p  <  0  then 
yields 


2  2 

<v>  p  /p  <1 
m 
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this  shows  that  <v>  decays  rapidly  for  p  >  p  ,  On  the  other  hand  one 

tn 

finds  that  p,..  is  always  about  3/4  p  ,  This,  together  with  Eq.  (14) 
U;  tn 

and  the  conclusions  derived  above  for  the  range  1.05  <  p  <  p ^  in¬ 
dicate  that  it  is  safe  to  assume  that 


0.60  $  <v>  s  1 


(31) 


as  long  as  p  £  p  . 

Hi 

In  the  limit  t  -*•  •  the  increase  of  n^  with  p  is  not  as  fast  as 

in  the  t  %  0  case  [compare  Eqs.  (3)  and  (7)].  On  the  other  hand, 

ip  has  a  slower  decay  because  S(t  •)  {  6(t  %  0),  so  that  n  remains 

e 

negligible  over  a  larger  p-range.  At  p  *  p(3)»  1;  although 

v(p^)  may  be  as  low  as  0.58  (because  nA  <  1)  we  have  p^j  >  p(i)» 
that  is,  p/0.  is  closer  to  p  .  For  the  lowest  p  values,  we  still 

\  ■>)  ra 

have  v(p  ■  1.05)  ■  0.60.  It  appears  that  the  inequalities  in  (31) 

are  still  valid  as  long  as  p  $  p  . 

m 

(c)  Bounds  for  arbitrary  t. 

According  to  Eqs.  (3)  and  (7),  we  have  n^(T  0)  >  n^(t  -*>  «*)  ■ 
There  is  therefore  a  net  decrease  in  the  ion  density  in  the  neighbor¬ 
hood  of  the  probe.  Let  us  assume  that  this  ion  depletion  is  monotonlc 
in  time.  Then 


6(t  %  0)  i  6(t)  i  6(t  -*  •) 


(32) 


since  6  goes  down  with  the  ion  density.  It  also  follows  that  for  any  t 
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0.60  <  <v>  <  1  ,  for  p  $  p_  .  (33) 

We  should  point  out  that  Inequalities  (32)  and  (33)  do  not  actually 
require  that  n^  be  monotonlcally  decreasing  In  time  but  that  the 
weaker  condition 

n±(T  %  0)  i  nt(T)  *  n^T  •*>  «)  (34) 

be  satisfied.  In  fact,  since  both  6  and  <v>  depend  on  n^  In  a  global 
way  only,  condition  (34)  could  be  locally  violated  without  Invalidat¬ 
ing  (32)  or  (33).  The  assumption  <  0  appears  well  justified 

at  least  In  a  global  way  since  the  current  density  j  exceeds  quite 
soon  and  remains  that  way  thereafter. 

We  now  define 


1/6  £  (1  +  0.63e2o?1./inJ;1/2)/6^1)  +  0.12 
(1)  P  z 

and  Introduce 


(35a) 


Y(c,  <|/p)  =  1/6  -  In  e 


(35b) 


Y  Is  given  graphically  In  Fig.  3.  We  shall  use  6  In  place  of  5(t) 
In  Eq.  (13).  The  maximum  possible  error  in  equating  6(t)  to  6 
goes  down  with  Increasing  and  e  We  find  that  the  error  Is 
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-1  -2 
less  chan  5%  for  c  ,  4>  ■  10  and  less  than  1Z  for  e  ■  10  ,  ip  z  10. 

P  P 

This  estimate  is  conservative  because  the  lower  bound  for  6(t)  has 

been  obtained  from  the  limit  t  ■+  “  while  only  a  limited  T-range  will 

be  considered.  Moreover  3iJ>/3p  will  only  be  used  under  Integral 

signs,  leading  to  a  smaller  error. 

We  shall  also  use  a  constant  value  v  instead  of  <v>  in  Eq. 

(13).  It  follows  from  (33)  that  putting  v  ■  0.80  produces  a  maximum 

possible  error  of  +  0.20.  It  should  be  noted  that  <v>  will  only 

appear  inside  Integrals  between  p  ■  1  and  p  i  o  :  this  will  obvlouslv 

m 

result  in  a  more  typical  error  for  <v>  %  v,  oi  +  0.10.  The  errors 
that  the  approximations  for  6  and  <v>  produce  in  the  determination 
of  the  ion  temperature  will  be  discussed  in  Sec.  VI. 

The  inequality  p  <  p^  will  appear  in  other  parts  of  our  analysis 
as  a  condition  for  the  accuracy  of  the  theory.  It  is  not  a  very 
precise  condition  in  the  sense  that  values  of  P  slightly  larger  than 
om  may  be  considered,  since,  as  just  indicated,  integration  over 

3tji/3p  reduces  the  error  in  <v>  a>  v.  o  may  be  considered  as  a  charac- 

m 

terlstlc  shielding  distance  since  it  is  close  to  which  is  itself 

very  close  to  exp  5  \  and,  moreover,  for  p  >  Pm  Che  plasma  becomes 
rapidly  quasineutral  and  both  potential  and  electric  field  become 
small.  It  may  be  illustrative  to  note  that  assuming  an  idealized 
charge  density  v  ■  1  up  to  a  certain  distance  Pp,  and  v  «  0  for 
o  >  Op,  the  first  two  members  of  Eq.  (15)  give  p_,  ■  p  .  The  first 

r  r  In 


and  last  members  then  give 
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6  -  (In  Pm  -  1/2) 
in 


-i 


which  can  be  used  as  an  approximate  equation  for  6, 


6  -  2(ln  <l>  +  In  6  -2  In  e  +  In  2  -  1) 
P 


-1 


(36) 


IV  -  The  overshoot  effect 

To  determine  the  region  A*  in  Fig.  2,  consider  the  radial 
motion  of  an  ion. 


2  2 
d_r  _  3V.  L 

^it2'  ±%9r  m.r3 


(37) 


The  angular  momentum  is  conserved  so  that  if  rQ,  are  the  Initial 
coordinates  of  the  ion,  L  ■  ro^vf^g  sin  In  non-dimensional  form, 
Eq.  (37)  reads 


c2  if  -  |*  +  26*  sin2  *  -f 
dr  p 


(38) 


2  - 1 

where  8*  ■  m,  (v*)  /2Z  <T  ;  from  Sec.  Ill,  3i(i/3p  ■  -  i|>  6o  (1 
l  *  i  e  p 

•  p  >.  Integrating  (38)  we  get 


-  2 
vp 


\  (e  i?)2  “  -  V  ln  p"  +  S*  (1  -  sin2  -|) 

0  p 


(39) 


2  v  ,  2  2. 

-  E  4  (P0  -  P  )  • 


1 
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For  (Oq,  it  )  to  belong  to  A*,  the  right-hand  side  of  Eq.  (39)  must 
be  non-negative  at  p  ■  1, 

S*o2  sin2  <j>0  s  8*  +  i|»p  5  In  pQ  -  e2v(p2  -  l)/4 
This  leads  to 

| sin  «J»0|  $  G(p0)/pq  ,  (40a) 

G  =  [1  4-  In  p2/a*  -  v(pj  -  l)/a*p2l1/2  (40b) 

where  a*  =  28*/i|>  6. 

P 

It  is  interesting  to  note  how  A*  changes  with  a*.  For  a*  >  1 

we  have  G(pq)/Pq  $  1  for  all  pQ  >,  1.  For  a*  <  1,  however,  there  is 

a  value  p  >  1  such  that  G(p  )/p  •  1;  then  G(p_)/p~  >  1  for 
q  q  q  ()  0 

pQ  <  p^.  As  shown  below,  the  regime  of  Interest  is  a*  <<  1;  con¬ 
dition  (c)  (Sec.  II)  ensures  that  this  is  satisfied  here.  The  upper 
boundary  of  A*,  P  ($) ,  is  then  given  by 

A 

(sin  $|  ■  min{l,  g<Pa)/pa}  »  (41) 

and  is  schematically  represented  in  Fig.  2.  As  a*  increases  point 
q  moves  down,  reaching  p  ■  1  at  a*  ■  1. 

The  preceding  analysis  shows  that  the  region  A*  is  determined 
by  considerations  of  angular  momentum  at  x  ■  0.  As  pointed  out  in 
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f 


Ref.  1  the  lower  the  Ion  temperature,  the  farther  from  the  probe  one 
must  go  to  find  ions  at  t  ■  0  that  will  miss  the  probe  due  to  their 
angular  momentum;  in  terms  of  Fig.  2,  when  8*  decreases  point  q  and 
the  whole  boundary  of  A*  moves  upward.  We  will  now  show  that 
J*/J  may  exhibit  a  significant  overshoot;  a  qualitative,  but  de¬ 
tailed,  picture  of  this  phenomenon  can  be  given  on  the  basis  of 
Fig.  2.  The  overshoot  may  be  seen  as  caused  by  the  sudden  set-up 
of  the  potential  field  which  traps  low  angular  momentum  ions  in  the 
neighborhood  of  the  probe;  in  the  steady  state  corresponding  to  a 
probe  potential  most  of  that  neighborhood  would  be  populated  by 
high  angular  momentum  ions. 

Up  to  a  certain  distance  from  the  probe  (say  p  <  r  / 3) , 

m 

3ip/ 3p  'v-  p  1.  For  Pq  in  that  region,  the  time  of  flight  to  the  probe 

T0(P0,  <f>0)  grows  (roughly)  linearly  with  p^  so  that,  according  to 
2 

Fig.  2 ,  a*  *v>  t  (except  for  t  <<  1)  as  long  as  p^  <  p^.  More 
precisely,  one  finds 

„  2  2  „  ,  7  2,  2 
a*  p  t  'v  ip  ot  /e 
m  p 


and  therefore 


Wj  ^  e  a*/Ti|>*^  'v-  $t/e 
"  P  P 


(42) 


j*/j  increases  with  both  ip  and  e  \  For  t  >  t  =  rrt(prt  ■  P  , 
m  p  q  U  u  q 

<Pq  *  tt/2),  da*/dt  will  be  nearly  a  constant,  because  a*  is  almost 


-24- 


llnear  in  p^;  thus  j*/jaB  will  taper  off.  For  even  larger  r,  particles 
from  outside  the  region  where  3<j*/3p  p  ^  will  begin  to  be  collected; 
the  field  now  decreases  rapidly  and  so  does  the  mean  velocity  of  the 
ion,  Pq/Tq.  Therefore  da*/dt  will  decrease,  and  j*/j—  will  peak  and 
begin  to  fall.  As  t  -*■  •,  j*/j  will  eventually  approach  unity.  In 
effect,  the  asymptotic  value  of  PA($)  sin  $  in  Fig.  2  is  obviously 

1/2 

p.  sin  $  -►  (*  /B*)  '  (43) 

A  p 

and  as  Pq  •  the  average  velocity  of  the  ion  becomes  v*,  so  that 


li"  Si*  .  2v?  r 

t  -  •  dt  x  pl  ,  ..2  1 


(44) 


.  lim  ..  . 

and  _  j*  *  j  • 

£-N»  J  J  00 

The  peak  in  j*/joo  is  sensitive  to  the  ion  temperature.  An  in¬ 
crease  of  £*,  and  therefore  a*,  produces  a  decrease  of  p^;  thus  the 
linear  growth  of  j*/j  given  in  (42)  is  Interrupted  at  a  shorter  time 
and  this  results  in  a  lower  peak  for  j*/.1  .  When  t*  reaches  unity, 

*  1  so  that  a*  never  grows  quadratically  in  r.  av  >  1  may  be 
roughly  considered  as  a  critical  condition  for  the  overshoot;  the 
numerical  results  of  Ref.  1  (largest  a*  %  0.70)  indicate  that  this  is 
actually  the  case.  For  decreasing  8*,  Increases  and  the  peak  of 
j*/j  goes  up.  F.ventually,  however,  a  value  of  8*  is  reached  for  which 
n  moves  out  of  the  3ijj/3p  •v#  p  *  region;  the  growth  of  1*/1  is  then 

q  oo 
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stop  ped  by  the  rapidly  decreasing  mean  velocity  of  an  ion  along  its 
trajectory,  and  not  by  the  fact  that  the  boundary  of  A*  is  reached. 
The  effect  is  the  same  for  all  smaller  B*  so  that  the  current  is  now 
Insensitive  to  the  actual  value  of  B*.  Since  Pm/3  grows  with  e 
it  is  clear  that  using  increasingly  smaller  probes  results  in  the 
possibility  of  measuring  increasingly  smaller  ion  temperatures.  We 
note  here  that  although  pm  grows  also  with  iff ,  so  does  p^  so  that  a 
change  in  the  probe  potential  has  no  effect  on  the  range  of  tem¬ 
peratures  that  can  be  measured  by  the  probe. 

To  determine  the  dependence  of  the  current  peak  on  the  various 
parameters  we  put  'v  p^/pm  in  (42);  according  to  Eq.  (41) 

p,  '  <£>1/2  «-  £V/2 


so  that 


f  *  <y2  ?n-  . 

This  shows  that  the  peak  is  fairly  Insensitive  to  the  value  of  e; 

this  prediction  is  supported  by  the  numerical  results  of  Ref.  1. 

The  current  becomes  insensitive  to  B*  when  p  'v  p  ,  or  t  'v  1;  this 

q  m  q 

roughly  corresponds  to  a  value  of  B* 

2 

B*  e 


(46) 
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The  maximum  peak  in  the  current  is  then 

J*/i  *  5/e  ,  (47) 

CO  p 

so  that  for  a  given  the  ratio  of  maximum  currents  for  two  different 

P 

values  of  e,  and  is  (Cj/e^Mln  tjVln  e^*) .  This  is  also  sup¬ 
ported  by  the  computations  of  Ref.  1. 

We  would  like  to  conclude  by  noting  that  (a)  these  qualitative 
results  can  be  obviously  extended  to  the  case  of  an  initial  ion 

Maxwellian  distribution;  (b)  for  o.  $  p  the  maximum  possible  error 

Am 

in  pa(4>)  due  to  the  error  in  putting  <v>  'v  v  =  0.80  is  around  1%, 
and  the  one  due  to  approximating  5  by  6  is  around  27.;  u)  the  inequality 
given  by  (40),  actually,  is  not  a  sufficient  condition  for  pQ  -  <(>0 
to  belong  to  A*:  For  5  ir/2,  the  right-hand  side  of  Eq.  (39)  must 

be  positive  for  all  p  between  p^  and  1;  in  other  words,  no  potential 
barriers  should  appear.  Defining 

x  -  o2/p2,  W  -  -  In  x  -  v  p2(l  -  x)/p2  +  a*(l  -  sin2  <Pn/x) 
u  u  m  u 

that  condition  reads 

-2  -2 

W(x)  >  0  for  Pq  <  x  <  1  whenever  W(p^  )  *  0.  (48) 

This  is  indeed  satisfied  when  p_  s  p  .  As  p.  increases  a  value  is 

Urn  U 

found,  dependent  on  both  a*  and  and  slightly  larger  than  p  ,  beyond 

u  m 

which  (48)  is  violated.  For  those  values,  however,  the  use  of  v  for 
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<v>  becomes  Invalid  so  Chat  no  clear  conclusion  can  be  reached  on  che 
matter  of  possible  transient  potential  barriers  for  (roughly)  pQ  >  pm. 
This  possibility  can  not  in  principle  be  ruled  out,  because,  for 
small  t,  the  somewhat  high  density  in  the  neighborhood  of  the  probe 
steepens  the  potential  field. 

For  ^  i  n/2,  the  sufficient  condition  is  that  W(x)  be  positive 
_2 

for  Pq  <  x  <  1  and  be  zero  at  some  x  >  1.  The  first  part  of  this 

condition  is  the  same  as  (48).  The  second  part  means  that  the  ion, 

which  for  >  it/ 2  moves  initially  away  from  the  probe,  must  turn 

back  at  some  p  >  p^;  this  condition  yields  a  new  boundary  for  A*, 

PA($).  which  is  also  represented  in  Fig.  2.  This  boundary,  however, 

plays  no  role  in  our  analysis  because  it  is  found  that 

Tn(P*($)»  >  t  .  where  r  is  the  largest  value  of  t  considered 

u  a  m  m 

(see  next  section  for  value  of  t  ). 

m 

V  -  Ion  current  computation 

The  derivation  of  a  formula  for  j/j  presents  no  essential  dlf- 
flculty  but  it  is  long  and  tedious.  We  will  now  give  schematically 
the  Mathematical  steps  Involved  in  the  derivation;  some  additional 
details  are  given  in  the  Appendix.  We  made  extensive  use  of  the 
assumptions  e  «  1,  ♦  »  1,  a*  «  1.  The  formula  derived,  however, 

is  not  just  asymptotic;  we  have  checked  throughout  the  numerical 
errors  of  the  various  approximations  or  simplifications  in  the  algebra. 
To  this  effect  we  specifically  assumed 

l(f 2  $  e  $  lO”1  ,  10  S  $  102  ,  a*  i  10_1 

P 
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(a)  Time  of  flight  to  the  probe  Tq(Pq,  $q)  . 

The  first  step  is  to  obtain  tq(Pq»  $q)»  the  time  of  flight  o> 
an  A*-ion  to  the  probe.  From  Eq.  (39)  we  get 

°°  p2 

'o  ■  7/2 1  dl>IB*(1  -  ;!  sl°2  V  -  v  ln  ^ 


-  o2>r1/2 


(49) 


for  <j>Q  *  n/2  and 


T0  "  To(n 


,)  +  2 


,1/2 


dp[8*(l  - 


0  .  2  .  V 

T8ln  V 


(50) 


for  ir/2  Hg  $  PT  is  the  zero  of  the  square  root  inside  the  in¬ 
tegral  and  represents  the  distance  at  which  the  ion  turns  back, 
heading  toward  the  probe. 

From  Eq.  (49)  one  obtains 


1/2  pQ  a*  Erf (In  pQ  +  a*/2)1/2  -  Erf(o£/2)1/2 

T0  “  9  P~  exp(F)  ~  ::  ’  2~7T  '  '  .  2  ~  .1/2 

m  (1  -  vgj^a*)  P0/om  -  o*  sin  <*>0 ] 


(51) 


where  a*  *  a*  cos  and 
c  0 


g.(a*)  -  3^2  exp  (a*) - - — ry= - 1. 


Erfc(3a*/2) 


1/2 


(52) 


Erfc(o*/2) 
c 
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From  Eq.  (50)  one  obtains 


To(ir 


$  )  +  -- 
* 0 '  o 


23/2  (a*)1'2 

- S.  u 


m  l  -  0*  -  vp_/o_ 
0  m 


1  -  2vp2/p2 

* - j  — j.  (53) 

3  C  (1  -  a*  -  »2/A2 
c  urn 


See  the  Appendix  for  details. 

(b)  Explicit  expression  for  Pq($q»  tq) 
Rearranging  Eq.  (51)  we  get  (for  $  ir/2) 


2  Vo  exp(-  °c)(1  '  °*  sin2  V 


0  ir{Erf (In  pn  +  a*/2)1/2  -  Erf (a*/2)1/2)2 
u  c  c 


where 


-  .  Erfc(ln  p  +  a*/2) 

•  U  +  *  T28.(a*)[l  -  — - — 


1/2 


-]"2)"1  (54) 


Erfc(<x*/2) 
c 


g2(a*)  *  exp  aV  ^Erfc  ^*c/2)1/2}  2 


(55) 


The  function  g2  remains  remarkably  constant, 

b 

0.732  -  g2(0)  $  g2(a*)  $  g2(0.1)  -  0.782  . 

2 

We  simplify  Eq.  (54)  by  approximating  exp(-  a*)  (1  -  a*  sin  <t>Q)  by 
1  -  a*,  g2  by  a  constant  g2>  and  the  bracket  by  unity.  The  first 
approximation  has  an  error  of  less  than  1%.  The  effect  of  the  second 
one  in  the  evaluation  of  the  ion  temperature  will  be  later  discussed; 
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we  choose  ■  0.75.  The  third  approximation  has  a  noticeable  error 
for  on  ■  0(1)  only;  because  of  the  large  value  of  p  .  pn  *  0(1) 

u  m  0 

corresponds  to  Tq  small,  in  which  case  the  term  with  the  bracket  is 
entirely  negligible.  We  get 


PplErf (In  pQ  +  a*/2)1/2 


Krf(a*/2)1/2}2  -  p2i2(l  -  a*)(ir  +  g0  v  t2)"1 


m  O' 


»2  *  ‘o' 

(56) 


Rearranging  Eq.  (53),  and  using  the  same  approximations  leading 
to  Eq.  (56),  we  get  (for  ir/2  $  $  $  ir) 


p2[Erf(ln  pQ  +  a*/2)1/2  +  Erf (a*/2)1/2]2  -  p2T2(l  -  a*)  (it  +  %2  v  t2)"1 


•  (1  -  A) 


-1 


(57) 


where 


a*  1/2  ®2  V  T0  -1 
A  -  4  Erf(“)A//  (1  +  U'  1 

C  T. 


)  *  [Erf (In  pQ  +  ~)1/2  +  Erf (  ~)1/2f2 


.  [—  Erf (in  pQ  +  ^-)1/2  -  r{l  +  Erf(^)1/2}  -  r2  Erf(^)1/2} 


U  g2  T0 


and  f  is  given  by 


[1  -  a*  -  v  g.(a*)  p^/o2!1^2 

9  2  ~  l 1 

1  -  a*  -  v  p  /p 
0  m 


,  a*(l  -  2v  p2Jp2)  a* 

_£  _C_  _  _  U  .H  __  C 

3  '  *  -  2  i  2.2  "  6 

(1  -  a*  -  v  p_/p  > 
u  m 


]  -  1 
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A  is  found  to  be  at  most  about  0.03  and  will  be  now  neglected.  This 
will  affect  the  value  of  a*,  to  be  computed  below,  in  little  more  of 
IX  since  A  appears  only  in  half  of  the  <t>^  range. 

Equations  (56)  and  (57)  may  be  now  substituted  by  a  single 
equation  valid  for  all 


r 2a*. 1/2 


p0u  -  (~-r'‘  cos  «j>0ni  - 


Erfc(ln  p„  +  a*/2) 
0  c 


1/2 


-] 


p  x_(l  -  a*) 
_m _ LL _ 


1/2 


1  -  (2a*/ir)1^2  cos  <f>0  (ir  +  v  g2  t^)1^2 


(58) 

1/2  1/2 

we  have  expanded  Erfc(a*/2)  for  small  (a*/ 2)  with  an  error  of 
less  than  1%. 

We  now  approximate  the  second  bracket  in  Eq.  (58)  by 

1  -  (2a*/ir)1^2  cos  <t>.  . /9 

1  -  Erfcfln  <o0  — ——ffi - 5,1  • 

(1  -  a*) 


This  involves  no  error  for  either  a*  ■  0  or  •>  ».  In  general, 
the  error  is  negligible  when  a*  is  small,  except  for  p^  about  one, 
that  is, tq  very  small;  then  our  analysis  fails  anyway  because  Eq. 
(5)  is  invalid.  Moreover  for  any  value  of  p^  the  error  vanishes  at 
some  4>q  between  0  and  n,  and  has  opposite  signs  for  larger  and 
smaller  angles;  this  results  in  a  partial  error  balance  in  the 
expression  for  a*.  Writing  then 

PQ[1  -  (2a*/ir)^2  cos  4>]  (1  -  a*)  *^2  ■  o 
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we  get 

2  _  o2(l  -  a*) _ 

0  (1  -  (2a*/n)1/2  cos  $]2 


(59) 


_  .1/2  .  .  -  -  2. -1/2 

a  Erf(ln  a)  ■  pm  i(n  +  v  ) 


(60) 


We  have  dropped  the  subscripts  In  both  $  and  t;  at  a  given  time  t, 

0q (4> ,  t)  represents  now  the  maximum  possible  p0  "collected"  for  each 
value  of  <J>.  Equation  (59)  exhibits  an  explicit,  remarkably  simple 
dependence  on  a*  and  $  (the  two  variables  over  which  further  Integra¬ 
tions  must  be  carried  out).  We  note  that  if  a*  ■  0,  p  equals  p  at 

IQ 

t  ■  2.80  for  all  If  a*  ■  0.1,  p  equals  pm  at  T  ■  2.80  for  <j>  »  76.4° 
and  at  smaller  (larger)  t  for  smaller  (larger)  4>.  Since  the  analysis 
Is  valid  up  to  distances  slightly  larger  than  pn>  we  may  roughly  state 
that  our  theory  begins  to  break  down  at  t  ■  Tm  %  2.80. 

(c)  Nondimensional  area  a*(i) 

To  obtain  a*(r)  we  integrate  over  p  up  to  Pq($»  t)  or  pa($)  which¬ 
ever  the  lowest,  and  then  Integrate  over  $.  In  Fig.  4  we  again  drew 
schematically  the  boundaries  of  region  A*;  we  also  included  several 
curves  of  the  p  ■  Pq($»  o)  family.  Notice  that  there  is  a  value 
such  that  if  a  <  o^»  p  ■  Pq($»  o)  and  p.  ■  p^($)  do  not  intersect.  For 
a  >  aL  there  are  two  Intersecting  points,  at  angles  $^(o)  and  $2(0); 

<J>1  *  $2  at  0  ■  oL»  From  Eqs.  (41)  and  (59)  we  get  (within  an  error 


of  less  than  IX) 
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C°s  * 


,2.1/2  in  &  -  1  1  +  a*  .1 

V  2  ~  1/2  - 

o  (a*) 


i  2 /  2  1/2 

In  o  /a 


(a*  -  h) 


1/2 


(61) 


,  rln  o2  2  ,ln  o2  -  1,1/2 
h  -  [— —  -  -  (  -  ) 

a  o 


<i  -  ^r1 

o  o 

m 


(62) 


obviously  is  given  by  a*  ■  h(o^) 

Let  us  denote  *  1(0^).  Then  for  t  <  we  get 


a* 


■h) 


kU)  -  *t» 


pQ(<i>,o) 


d<t> 

J  J 

0  0 


p  do 


Pq (d> »  0)  dd> 


(63) 


tj/2  + 


cos 


IT 

r 


d$ 


2 

Pn($  » 


P) 


The  area  a5^  corresponds  to  the  probe  cross  section  and  to  that 
part  of  the  shaded  region  lying  inside  the  curve  p  *  Pq(<{>,  k°th 
were  included  in  for  simplicity.  The  point  p^,  is  the  inter¬ 

section  of  the  line  p  sin  $  ■  1  (u  $  *  n/2)  and  (59).  We  find 


a2  (1  -  a*)  dd» _ _ 

[1  -  (2a*/ir)1^2  cos  p]2 


it  o2(l  -  a*) _ 

(1  -  2a*/it)3/2 


(64) 


the  dependence  of  Kq.  (64)  on  a*  amounts  to  less  than  1%  and  will  be 


neglected. 
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a(l)  “  na 


(65) 


Within  the  same  error  we  find,  after  a  long  transformation, 


.  2  .  -1  1  .  ,  2  1Nl/2  _  ,2a*xl/2  .... 

a(2)  "  2  +  °  sin  a  +  (a  “  1)  “  2(-~ )  a  .  (66) 


For  T  >  tl»  ^(2)  unchanged  hut  a5^  becomes 


^(a 


d>2  (cr ) 


^tl)  "  {  °0(*’  °)  d<>  +  j  (<!>)  d$  +  |  p2($,  a) 


♦  1(o) 


<t>2(o) 


The  first  and  third  integrals  give 


(67) 


oQ  d<Ji  + 


2  »A 
°0  d* 


_ C2a.Vn),1/2  (1-  a*)  o2  sin  f _ 

(1  -  2a*/ir)  [1  -  (2a*/n)1/2  cos  $1 


1/2 


.  (1  +  “lffl  +  (2a*/w)  1 1/ 2  i, 

+  “  .  .  3/2  (2  +  tn  l{  T/2}  tn  2 


(1  -  2a*/ir) 


1  -  (2a*/ir) 


1/2 


and  after  substitution  of  $^(a)  and  ^(p)  and  some  tedious  transforma¬ 


tions  we  obtain 
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|  o2Q  d*  +  j 


2  j,  2r _  ,  -1  ,h  >1/2  .  ,  .  .,1/2  ,1/2 

0Q  d<j»  ■  a  [2  sin  (^)  +  (a*  -  h)  h 


^  4  + _ 12  8  vo 


*  ,2  2.2  .  2,2 

it  In  a  ir  p  In  a  tr(ln  o  ) 


})  . 


(68) 


correct  up  to  1Z.  When  h  -*•  a*,  Eq.  (68)  approaches  ira  and  d>  -►  <(>  ; 
thus  we  recover  Eq.  (65). 


To  calculate 


p*  d4>  we  first  obtain  an  explicit,  approximate 


.  *1 


expression  for  (<p)  from  Eq.  (41) 


p2  *  (a*  sin2  $  +  v/p2)  *  Inf (a*  +  v/p2) 


m 


(a*  sin2  cp  +  v/p2)  *  1.32  p2  exp  a*] 


(69) 


where  p  =  p.(ir/2)  is  given  by 
q  A 


(a*  +  v/p2)  p2  31  In  p2  +  a*  +  v/p2 

*n  q  q  m 


The  maximum  error  in  Eq.  (69)  is  found  to  occur  at  the  largest  a*. 

and  amounts  to  7.7%.  The  following  points  should  be  noted  however: 

(1)  That  error  occurs  at  a  certain  value  of  $  so  that 

*2 

—2 

p.  d<J>  will  obviously  be  more  accurate  f the  more  so  since  the  error 

i  A 
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in  Eq.  (69)  can  have  both  signs].  (2)  After  the  integration  is  per¬ 
formed  we  reverse  the  change  from  Eq.  (41)  to  Eq.  (69),  substantially 
reducing  the  error.  (3)  The  error  is  maximum  when  the  contribution 


oA  d<j>  to  the  final  expression  for  the  current  is  minimum.  Carry¬ 


ing  out  the  integration,  and  substituting  and  ^  from  Eq.  (61) 
we  finally  obtain 


-2  2  r  „  ,  -1  ,h  ,1/2,,/.  .  v  ,-1  .  v  , 

°a  «* ["  ■ 2sln  <■„*>  lia  +  — 2>  +  2 

a*  o  3a*  o 

m  m 

-  ~  [h(a*  -  h)]1/2  U  -  -2  (-1-n-0-2-^)2} 

In  a* 


+  2 (  -  X>1/2  <a*  “  JO 1/2  [In  o2  (2  -  ~)  +  1 

.  4  TT  Z 

In  o 


2  .  4  .  C  2  —  In  a  -  2/3, 

-  7 +  '  2  +  i*  -  3  v  -2—2 1  •  (70) 

tt  In  a  r_(a*) 

m 

«„  »2  -  2)  [1  +  f  (1  -  -i-J  +  \  -4  -V--- 

In  a  (In  a  )  p„  In  o 

m 


(d)  Nondimensional  time-averaged  current  j/jo 
From  Eqs.  (8),  (10)  and  (11),  it  follows  that 


=  t  iil'I  dr' _ -E— ; - f  da- 

T  i  3-  (8*  )1/2  T  l  a 


exp(-  j— )  a*  ; 


? 
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(notice  that  m^(v*)  /2ieT^  ■  3*/6  ■  a*/a) .  The  integral  in  Eq.  (72) 
may  be  rewritten  as 


|  da*  exp(-  a*/a)  a* 
0 


h(0) 


|  da*  exp(-  a* /a)  a^^ 


|  da*  exp(-  a*/a)  a*^  -  j  da*  exp(-  a*/a)  a*2^;(73) 


h(cr) 


0 


Eqs.  (65)  and  (67)  must  be  used  for  the  first  and  second  integral 
respectively.  Notice  that  values  of  a*  much  larger  than  0.1  are 
considered  in  Eq.  (72);  if  a  $  0.1,  however,  our  estimate  of  the 
errors  produced  by  those  simplifications  based  on  the  assumption 
a*  $  0.1  still  holds,  because  as  a*  increases  beyond  0.1,  the  in¬ 
crease  in  those  errors  is  much  slower  than  the  exponential  decrease 
of  the  exp(-  a*/a)  factor.  For  the  last  integral  in  Eq.  (73)  we 
find 


|  da*  exp(-  a*/a)  a*2^  ■  a[u/2  +  o2  sin  1(l/o)  +  (o2  -  l)1^2  -  (2a)*^2  a]. 


For  the  first  and  second  integrals  we  find  (after  considerable 
simplification) 


(74) 


da*  exp(-  c.*/a)  “  auo2[(l  +  2n  *)  Erf  n  ^2  -  2n  *  +  2 (tt n )  *^2 

•  exp(-  n  *)  -  (In  a2)  *  (4(itn)  1^2  exp(-n  *)  Erfc  n  *^2 

2  1/2 
_  -1  „  -1.,  2  v  a*  ,  . -1/2  ,  -1>  ,  1  +  4/tt  n  ,  -1N 

2n  t1(2n  ))  -  j  ~~2 — ~~2  exP(“  n  )  +  ~\J2 - 2~  exp^‘n  ) 

p  In  a  no 

m 


(75) 
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correct  up  to  1%.  We  have  deflnfid 


a  a 


In  ad 


(76) 


Our  results  can  now  be 


summarized  as  follows: 


j-  (t-  E'  s.  O  -  ~’~c  ?*  — 


p’  (sy V*  7  X(°’  "•  1 


on 


x  ■  X1  *  X2  +  x3 


(78) 


X1  *  Cl  +  2n_1)  Erf  _  ,_-l 


-  2n  *  +  2(nn)~1/2  exp(-  n-lj 


1  - 


-2 


r„>  t4("p,'1/2  «*-  n-i>  wc  ,-»/* .  2n-i  V2n-i)}_ 


X  -  IL.{2  2  V  02  -  1  _i/o 

ffl/2  a2  ’  3  n2  ln  2  (,rn)  }  exp  (-  r f1)  , 


n_  In  0 
ID 


X3  ■  -  -  r  sin'1  -1  -(j L--  °~V1/2 

2o 


3  -  2  ~  *  si"  ;•  - - (i  -  (lUn/jl^j 


iro 


”here  n*  “m  are  8lven  i„  terms  ot  t.  e,  B 


o2  -  1 


and  \i)  by 
P 


(79) 


(80) 


(81) 


°m  '  (2*n  «1/2/c  , 


(82) 


1/2 


-  CmT(ff  +  V  g  T2)'1/2  , 


O  Erf (In  a) 


(83) 
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n  -  26(o2  -  l)/^i  6  In  a2  (84) 

P 

and 

v  -  0.80,  I-  -  0.75,  S  _1-  In  e_1  +  Y(e,  ip  )  (85) 

*•  P 

where  Y  is  given  by  Fig.  3. 

The  instantaneous  nondimens ional  current  j/j  can  be  obtained 
from  Eqs.  (77)  -  (85)  and  identity 


~  (T  j~) 


Our  emphasis  on  j/jm  is  based  on  (1)  the  greater  simplicity  of  the 

analysis  of  a*(t)  as  compared  to  that  of  da*/dt,  and  (2)  the 

possibility  of  using  j/j  ,  but  not  j/j  ,  in  an  alternative  method 

of  ion  temperature  determination  (see  Sec.  VI). 

For  comparison  with  the  approximate  numerical  computations  of 

-2 

Ref.  1,  we  give  in  Fig.  5  j/j  versus  t  for  e  *  10  ,  ip  *  15  and 

00  p 

several  values  of  8,  from  both  those  computations  and  F.qs.  (77)  - 
(85).  The  agreement  is  good. 

o 

X  is  graphically  given  in  Fig.  6  for  5  <  o  <  5.78  x  10"  (this 

-1  -2 

is  the  largest  value  of  o  in  the  range  t  5  2.80,  10  >  e  >  10  , 

2 

10  <  \|»  <  10  ) .  We  notice  that  X  approaches  X(o  -*■  ®)  quite  rapidlv: 

is  never  more  than  2%  of  X,  X^  decays  quite  rapidly  as  o  increases, 
and  the  last  term  in  Eq.  (79)  is  always  a  small  fraction  of  X^ .  We 
also  note  that  Fig.  6  does  not  depend  on  cn'  we  used  an  appropriate. 


I 
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intermediate  value  of  o  for  each  curve  in  that  figure,  and  the  result- 

tn 

ing  error  for  other  values  of  c  is  less  than  +  1%. 

m 

Finally  we  point  out  that  to  present  our  final  forrulae,  Eqs. 

(77)  -  (85),  we  substituted  In  o^/(l  -  a  ^)  for  In  in  Fqs.  (75)  and 

(76);  otherwise  j/j  would  go  to  infinity  as  a  -*■  1  (t  -*>  0) .  Any 

-2 

function  F(o)  could  be  used  instead  of  (1  -  o  )  as  long  as 

2  -2 
In  o  /F(o)  *  0(1)  at  o  =  1  (the  actual  choice  F  *  (1  -  a  )  originates 

from  an  analysis  of  <f> ^ (o )  and  ^(o)  at  a*  ■  1].  This  arbitrariness 
justifies  the  procedure,  and  is  due  to  the  smallness  of  the  contribu¬ 
tion  of  the  a*/o  >>  1  range  to  the  Integral  in  Eq.  (72).  We  also 
2  1/2 

wrote  (o  -  1)  for  o  in  the  last  term  of  Eq.  (74);  otherwise  we 

again  would  have  j/jw  +  "  as  t  >  0,  for  a  i  0.  The  error  is  due  to 

the  simplification  leading  from  Eq.  (58)  to  (59).  We  stress  that 

these  corrections  are  only  required  for  x  0;  this  regime  was 

explicitly  excluded  from  our  analysis  because  of  the  failure  of  Eq. 

—  1/2 

(5).  Still  note  that  with  those  corrections  j/j^  (B/2^p)  as 

1/2 

x  ■*  0,  while  the  true  value  is  (ir/2)  larger,  an  error  of  only 

20%. 

VI  -  Conclusions 

Analytical  formulae  have  been  derived  for  the  transient  ion 
current  to  a  long,  cylindrical  Langmuir  probe  in  a  quiescent, 
collisionless  plasma,  after  a  switch  in  the  probe  potential  from 
zero  to  a  large  negative  value  (-  eVp  >>  icTe,  tcT^);  the  initial  ion 


-41- 


distribution  function  is  assumed  to  be  Maxwellian,  and  the  electron 

Debye  length  1..  to  be  large  compared  with  the  probe  radius  r  . 

D  p 

The  transient  current  density  j(t)  exhibits  a  large  overshoot  be¬ 
l/2  -1 

fore  approaching  its  new  steady-state  value,  1*2  n  e  No 

00 

1/2 

•  (-  Z.  eV  / m.)  .  The  overshoot  is  sensitive  to  the  ion  tem- 

i  p  i 

perature  and  may  be  used  for  diagnostic  purposes.  Assume 

-  eV  /icT  =  $  and  r  /A_  =  e  to  be  known  (from  conventional  steady 
p  e  p  p  D 

state  measurements,  say)  aud  let  j/j  be  the  experimental,  non- 
dimensional,  time-averaged  current  density  determined  as  a  function 
of  t  =  Upt,  where  Up  is  the  ion  plasma  frequency.  A  single  point  in 
the  curve  j/jm  versus  t  can  yield  the  ion  temperature:  For  any 
couple  of  values  (j/jM,  t),  the  quantities  o  and  X,  as  defined  by 

o  Erf (In  o)1/2  -  (2*p  5)1/2  e'1  T (n  +  v  g 2  x2)“1/2  ,  (86) 

X  -  (j/jj  (8^p)1/2  (^e)"1  (x/o2)  ,  (87) 


may  be  evaluated;  the  quantity  n(o,  X)  may  then  be  obtained  from  Fig. 
6,  and,  finally,  the  temperature  ratio  from 


ZiTe 


1 1>  6  In  a 


(88) 


In  these  equations 
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v  -  0.80  ,  g2  -  0.75  , 


(89) 


6  1  -  In  e_1  +  Y(e,  <p  ) 

P 


(90) 


and  Y  is  given  by  Fig.  3.  It  should  be  noted  that  can  be  determined 
even  if  Tfi  is  unknown.  In  effect,  we  can  write  Eq.  (88)  as 


_  „  „  6  In  a 

kT.  -  -  Z.eV  — - -  n 

1  1  P  o2  2(1  -  o'2) 


1/2  —  2  — 
ili  /e,  and  therefore  6/a  ,  is  not  a  function  of  T  *5  does  depend 
P  c 

on  the  electron  temperature  but  very  weakly:  the  logarithmic  de¬ 
pendence  in  In  e  1  is  partially  balanced  by  Y.  Thus 


Ti  *  ln-  n[o,  X<°>1 

l-o 


a  *  [<5(Te)] 


1/] 


in  the  range  of  interest  of  the  theory  [large  o,  1  <  n  <  100],  no 
sensible  dependence  of  on  Tg  is  found. 

/ 

In  the  process  of  deriving  formulae  (86)  -  (90),  certain 
functions,  6,  <v>  and  g2>  were  found  to  vary  between  narrow  bounds; 
they  were  approximated  inside  a  number  of  integrals  by  constants 
6,  v  and  g2  as  given  by  Eqs.  (89)  and  (90).  It  is  now  possible  to 
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discuss  the  errors  in  that  may  result  from  these  approximations. 
For  large  a  we  can  write 


Tj.  *  (ir  +  g2  vt2)  n(X  (a)  ] 


a  *  61/2  (it 


.  -  -  2. -1/2 
+  g2  ) 


and  for  2  <  n  <  30  we  roughly  have 


X  %  1  -  0.20  In  n 


From  this  and  the  bounds  found  earlier  for  6  we  get  an  estimate  of 
the  maximum  possible  error  due  to  the  indeterminacy  in  6  of  less 
than  10%;  when  t  =  1.5  a  similar  result  is  found  for  the  product  g2v. 
We  would  like  to  stress  that  these  are  very  conservative  criterions 
and  that  the  actual  error  should  not  be  larger  than  a  few  percent. 

We  note,  however,  that  the  possible  error  due  to  g^v  grows  with  x 
and  reaches  20%  at  x  ■  2.80.  Again  this  is  a  very  conservative 
estimate,  but  clearly  indicates  that  the  theory  should  not  be  used  for 
x  larger  than,  say,  3. 

There  are  certain  experimental  conditions  that  must  be  met  if 

this  method  of  measuring  ion  temperature  is  to  be  useful.  In  terms 

of  the  four  independent  parameters  x,  e,  B  =  T./Z.T  and  <i>  governing 

lie  p 

the  function  j/j  ,  the  validity  of  the  present  theory  is  roughly 
restricted  to  the  range 
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0.1  <  t  <  3  ,  10"2  $  e  i  10-1 

^  <$(c,  ij)  )  . 

3  5  -E— 20 - *  10  5  *  10  * 

-2  2 

Conditions  e  >,  10  and  i|*  <  10  could  readily  be  relaxed  in  the 

P 

theory.  Conditions  iji  >10  and  20B  $  iji  5  could  be  somewhat  relaxed, 

P  P 

but  any  substantial  relaxation  would  complicate  the  analysis 
extraordinarily;  these  requirements,  moreover,  are  easily  met 
experimentally . 

Condition  e  <  10  ^  or,  more  generally,  that  e  be  small,  is  a 
more  basic  one,  since  when  e  increases  beyond  10  \  say,  the  ion 
current  becomes  rapidly  insensitive  to  and  the  method  ceases  to 
be  useful.  More  specifically,  note  that  for  n  <  1,  j/j  becomes 
independent  of  B  (see  Fig.  6)  so  that  from 


2 

T _ 

TT  +  0.6  T2 


it  appears  that  for  a  given  e,  there  is  a  value  of  B  below  which  the 

method  can  not  measure  ion  temperature;  this  value  is  approximately 

-4  -2  -2  -1 

2  x  10  for  e  *  10  and  10  for  e  *  10  .  Since  a  probe  of  radius 


r  ■  10  cm  can  readily  be  obtained  by  conventional  methods,  re- 

-1  -2  5 

quirement  e  $  10  is  met,  at  least,  by  plasmas  with  z  10  cm. 

The  requirement  is  more  easily  satisfied  at  higher  electron  tem¬ 


peratures  and  lower  electron  densities 
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The  restricted  range  0.1  <  t  <  3  creates  another  experimental 

limitation.  To  resolve  time  on  the  ion  plasma  frequency  scale  it 

is  necessary  that  the  time  constant  KC  of  the  measuring  device  be 

very  small  compared  with  the  ion  plasma  period,  2 tt / .  Assuming 

that  the  resistance  R  is  governed  by  the  requirement  that  the  probe 

-3 

current  give  rise  to  a  potential  V  *  y  *  10  volts  across  it,  and 

-12 

that  the  capacitance  of  the  circuit  is  C  •  x  .  10  farads,  we 
obtain  the  condition 


-3  1/2  3 

xy  «  3  •  10  J  eU/XD)  (No  \*)  (91) 

where  l  is  the  length  of  the  probe.  For  typical  instrumentation,  y 
should  not  be  less  than  about  unity,  and  x  would  be  larger  than  unity, 
but  of  that  order.  Note  that  in  the  above  estimate  we  used  j  in 
computing  R  -  V/2irr  i,  j;  since  J  >>  j  ,  Eq.  (91)  is  conservative. 

p  00 

A  final  consideration  to  be  made  is  whether  the  displacement 
current  due  to  charge  accumulation  on  the  probe  will  mask  the  true 
current  (particle  flux)  which  is  to  be  measured.  Consider  the  charge 
balance  for  the  probe.  If  □  is  the  surface  charge  density,  we  have 


jD  =  3 %/n  -  j  -  le  -  ,iT 

where  j  and  j  are  the  ion  and  electron  particle  current  densities 
and  2Trrp&  is  the  current  measured  in  the  circuit  outside  the  plasma. 
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jg  may  be  neglected  because  it  is  exponentially  small  in  ^  so  that 
our  method  will  be  useful  as  long  as  is  small.  The  displace¬ 

ment  current  j  can  be  rewritten  as 


1_  3_  /“  _3V . 

4 it  3t  '  3r  'r-r 


The  function  6(t)  was  studied  in  Sec.  Ill;  its  variation,  and  there¬ 
fore  j  ,  is  due  to  the  rearrangement  in  space  charge  density  follow¬ 
ing  the  application  of  the  potential  to  the  probe.  Time  averaging 
we  find 


By  x  %  0  we  mean  a  short  time,  t  <<  1,  after  which  the  electrons  are 

in  equilibrium  with  the  field  and  6(x)  changes  very  little;  excluding 

this  time  (of  the  order  of  a  few  electron  plasma  periods)  in  the 

averaging  reduces  considerably  the  value  of  j_/j  . 

D  00 

A  conservative  estimate  of  [5(t)  -  6(t  %  0)]  t  ^  would  appear  to 
be  [6(»)  -  6(x  %  0)  ]  (2tt)  ^ ;  the  calculations  of  Ref.  1  indicate  that 
the  asymptotic  approach  to  the  final  steady-state  value  requires 


several  ion  plasma  periods.  With  the  above  estimate  we  find  that 


—  2  -2 

j^/j^  is  just  less  than  one  for  \p  *10  and  e  ■  10  ,  and  goes  down 

-1  —  2 
with  both  \p  and  e  .  Since  j/j  is  typically  of  order  of  10  for 

P  00 


2  -2 

ip  *  10  and  c  *  10  ,  and  varies  with  \p  and  e  in  a  way  similar  to 

P  P 
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! 

j, 

j_./j  ,  it  appears  that  j-./j  is  here  of  order  of  a  few  percent.  This 
validates  the  method. 

Even  if  j^/j  were  not  small,  the  method  could  be  essentially 
salvaged.  We  propose  the  application  to  the  probe  of  a  potential 
pulse  of  width  T  %  3/u^,  Instead  of  the  step-wise  potential  discussed 
until  now.  If  the  time  average,  now,  extends  from  a  few  electron 
plasma  periods  before  the  beginning  of  the  pulse  up  to  a  few  electron 
plasma  periods  after  the  pulse  is  ended,  it  is  easy  to  verify  that 
the  displacement  current  cancels  out. 
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Appendix 


To  carry  out  the  integration  in  (49)  we  rewrote  this  equation  as 


T  „  0-1/2  P0  f  ..  ds(l  -  P) 

0  p_  ]  _l/2 


-1/2 


m  '  0  s  (-  In  s  +  a*) 
-L  c 


1/2 


(Al) 


where  p 


P1  +  p2  and 


a*(l  -  s)  sin 

*1  s(-  In  s  +  a*)  * 

c 


v(l  -  s)  p 


0 


2  p2(-  In  s  +  a*) 

m  c 


From  Eq.  (Al)  we  immediately  get 


To"  11 


1/2  °0 


e*p(~)  [Erf  (In  pQ  +  ^-)1/2  -  Erf(^)1/2][1  -  p]’1/2  (A2) 
pm 


—  -2 
where  p  is  the  value  of  p  at  some  appropriate  s  between  and  1. 

To  obtain  p  we  expanded  the  bracket  in  (Al).  In  the  range  of  variables 

2  2  3 

considered  we  found  that  the  p^,  p^f  Pj  and  P2  terms  had  to  be  re¬ 

tained  to  achieve  an  error  of  about  1%.  Within  an  error  of  (at  most) 
a  few  percent,  all  the  resulting  terms  in  the  expansion  of  could  be 
summed  to  obtain  an  expression  like  (A2)  where 


p  *  a*  sin2  $0  +  v(p2/p2)  g7 (a*) 


g^  being  given  by  Eq.  (52).  The  result  was  Eq.  (51). 
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For  the  second  integral  in  Eq.  (SO),  the  expression  inside  the 
square  root  could  be  approximated  by  a  polynomial  of  second  degree. 
The  result  was  Eq.  (53). 
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By  t  %  0  we  mean  a  time  t  <<  1  such  that  Eq.  (5)  is  already 

satisfied,  while  the  ions  have  had  no  time  to  respond  to  the  field. 

2  2 

It  will  suffice  to  notice  that  i(i(t  -  ®,  p  >  P^)  >  PP^/p  »  as 
follows  from  the  condition  for  orbital  motion  limited  current 
collection. 

1/2  -4 

If  6  >.  0.1,  and  a  probe  with  rp  ■  10  x  10  cm  can  be  obtained, 

-3 

one  can  move  this  condition  down  to  *  10  cm. 


of  ions  that  will  be  collected  at  positive  times 


vs.  n  and  o.  Curves  end  at 


